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ABSTRACT 



The first Search for Extra- Terrestrial Intelligence (SETI) conducted with Very 
Long Baseline Interferometry (VLBI) is presented. By consideration of the basic 
principles of interferometry, we show that VLBI is efficient at discriminating 
between SETI signals and human generated radio frequency interference (RFI). 
The target for this study was the star Gliese 581, thought to have two planets 
within its habitable zone. On 2007 June 19, Gliese 581 was observed for 8 hours 
at 1230-1544 MHz with the Australian Long Baseline Array. The dataset was 
searched for signals appearing on all interferometer baselines above five times 
the noise limit. A total of 222 potential SETI signals were detected and by 
using automated data analysis techniques, were ruled out as originating from 
the Gliese 581 system. From our results we place an upper limit of 7 MWHz -1 
on the power output of any isotropic emitter located in the Gliese 581 system, 
within this frequency range. This study shows that VLBI is ideal for targeted 
SETI, including follow-up observations. The techniques presented are equally 
applicable to next-generation interferometers, such as the long baselines of the 
Square Kilometre Array (SKA). 

Subject headings: Techniques: interferometric, Radio continuum: planetary systems, 
Stars: individual(Gliese 581) 
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Introduction 



The Search for Extra- Terrestrial Intelligence (SETI) seeks evidence for life in the 
Universe, through the detection of observable signatures from technologies that are expected 
to be possess ed by advan ced civilizations. Such searches are mainly co nducted at radio 



wave lengths ((Drake 
2002 ) . 



20081 ) and, to a lesser extent, optical wavelengths (IReines fc Marcy 



The goal of SETI experiments conducted at radio wavelengths is to detect signals that 
are intentionally aimed at broadcasting the civilization's existence to othe rs, and/or sign als 



that unintentionally leak from the civilization's communications system ( ISiemion et al. 



2008 



Fridman 



20111 ). Intentional signals are exp ected to be narrow in frequency 



1 kHz), while unintentional signals may not be ( iTarter 



2001 



Fridman 



201lh . 



2004 ; 



in frequency ( 


1 Hz to 


Siemion et al. 


2008 





An ope n question is the frequency 



intelligence (ITarter 



2004 



Shostak 



2009 



range over whi ch to search for extra-terrestrial 



Fridman 



20111 ). The frequency range from 1 



to 10 GHz i s generally p referred, however after decades this range is still very much 



unexplored (Tarter 



( ITarter 



2001 



20041 ). For leakage signals, the lower end of this range is preferred 



20041 ). For deliberate signals from extra-terrestrial intelligence, the most 



common suggestions are the 1420 MHz hydrogen line and multip 



7T x 1420 MHz (4462.336275 MHz) (Blair et al. 



1992 



Harp et al- 



es of 1420 MHz, including 



2010 . 



SET 

searches 



experiments are divided into two categories: (1) sky surveys and (2) targeted 



Shostak 



2009j). The sky surveys make no assumptions about the location of 



the SETI signals and perform raster scans of the observable sky. Current sky surveys 



are SETOhome, ASTRO 



telescope in Puerto Rico (jWerthimer et al. 



'TJLSE and SER 



)NDI 



2001 



V conducted with the 300 m Arecibo 



Siemion et al. 



20081 ) . Targeted searches 



sequentially examine specifically chosen stars that are deemed to have a chance of 
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harbouring a planet that could sustain life (ITarterl I2004J ; lFridmanll201ll ). A list of 17,129 
stellar systems that are potentially habitable to complex life forms was compiled by 



Turnbull fc Tarter! (120031 ). 



The currently active Kepler Space Missio n aims to search for Earth-sized planets in 



and near the habitable zone of Sun-like stars ( Borucki et al. 



20101 ). As of February 2012, 



Kepler has found over 2,30 planetary candidates, with ~46 candidates within habitable 



zones (jBatalha et al 



20121 ) . As Kepler begins to confirm these planets and proceeds to 



discover more Earth-like planets in habitable zones, these planets will be investigated for 
extra-terrestrial intelligence. 

This paper explores the suitability of using Very Long Baseline Interferometry (VLBI) 
for targeted SETI from potentially habitable planets. SETI is listed under the The Cra dle 



of Lif e Key Science P roject for the Square Kilometre Array (SKA 



( 120041 ): 



Tarteri (120041 ): 



Schilizzi et al 



Carilli fc Rawlings 



(j2010f )). This paper also aims to set a foundation for 



future VLBI SETI projects, including the use of the long baselines of the SKA. 

In this paper we present the first SETI experiment conducted with VLBI. Section [2] 
introduces the target for this study, Gliese 581. Section [3] discusses the advantages of using 
VLBI for targeted SETI and describes the technique of using VLBI for SETI. Section @] 
describes the observations and data analysis, and presents the results of the experiment. 
Discussions and conclusions are presented in Sects. |5]and 6, respectively. 



2. Target 



The target for this pilot study is the M-dwa rf star Gliese 581 (G1581), located 20 



light-years (ly) distant in the constellation Libra (jUdry et al 



2007 



Vogt et al. 



20101 ). 



In 2007, 



Udry et al 



(120071 ) announced the discovery of a planet on the edge of G1581's 
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habitable zone. The planet, Gliese 581d (G1581d), is the fourth planet in the Gliese 581 
system. It has a mass of 5 Me, an orbital period of 83 days and an orbital semi-major axis 



of 0.25 AU (~40 milliarcseconds (mas 



■fl 



Wordsworth et al. 



( 120101 ) discussed the effects of 



surface gravity, surface albedo, cloud coverage and CO2 density on the surface temperature 
of G1581d. Their results showed that G1581d may have the necessary conditions to sustain 
liquid water on its surfac e. Thus, G1581d may be the first confirmed exoplanet with the 



possibility to sustain life (IWordsworth et al. 



20101 ). 



Vogt et al. 



(120101 ) claimed to have discovered a planet within the habitable zone 
of G1581. The suspected planet, Gliese 581g (G1581g), is located between G1581d and 
Gliese 581c (the third planet). G1581g has a proposed mass ~ 3 Me, an orbital period of 
37 days and an orbital semi- major axis of 0.15 AU (projected angular distance ~23 mas). 
The ex i stence of a low mass planet in the habitable zone of G1581 has been questioned by 



Tuomi 



(2011). 



No radio emiss ion was detected from the position of G1581 by the VLA FIRST survey 



( IBecker et al 



19951 ) above a limiting flux density of 1 mJy at 20 cm. 



3. Using VLBI for Targeted SETI 
3.1. VLBI 

Very Long Baseline Interferometry (VLBI) allows, by the combination (via a correlator) 
of signals from multiple radio telescopes, the emulation of a telescope the size of the 
maximum telescope separation, which is generally hundreds to thousands of kilometres. The 
outputs of the correlator are 4-dimensional (baseline, time, frequency, polarisation) datasets 

1 this is the projected angular distance as seen from Earth 
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consisting of multiple spectral channels and time averages over ~1 second. The primary 
characteristic of VLBI is its very high angular resolution, which falls in the milli-arcsecond 
regime, the highest resolution in a stronomy. This high re solution makes VLBI ideally suited 



for high precision astrometry (see 
energetic compact objects. 



Thompson et al. 



2001 



chap. 12), and observation of high 



Existing VLBI instruments are able to monitor a wide 



10 MHz (the Low F requency Array, L OFA R: 



Horizon Telescope: 



Broderick et al 



Stappers et al. 



range of frequencies, from 



(boilh ) to 230 GHz (The Event 



( 120111 )). In addition, current VLBI systems such as 
the LBA and the European VLBI network (EVN) include some of the most sensitive radio 
telescopes in the world. 



To date there has been little or no use o 



project SETI-Italia (jMontebugnoli et al. 



VLBI Medicina telescope. SETI-Korea (IRhee et al. 



VLBI in SETI investigations. The SETI 



200 lj) uses data commensally from the 32 m 



20101 ) has proposed the use of existin 



Slysh 



and n ew telescopes from the Korean VLBI Network (KVN) for SETI. In addition 
( 1199 it ) proposed to use VLBI or space- VLBI satellites to measure the angular size of 
suspected SETI radio sources. Pseudo-interferometric observations were conducted by the 
SETI Institute, using both the Arecibo Tel escope and the Lovell Telescope, Jodrell Ba nk 



Observatory between 1200 and 1750 MHz (iBackus fc The Project Phoenix Team 



2002|). 



The survey searched for signals from nearby stars with appropriate frequency drift and 
offset caused by the rotation of the Earth. 



3.2. Radio Frequency Interference (RFI) Mitigation 



Radio frequency interference (RFI) is human-made, unwanted radio emission within 
the telescope's field of view. Potential sources include mobile phones, radars, television 
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signals, FM-radio, satellites etc. Separating RFI from a real SETI signal is one of the most 
challenging problems faced by SETI projects. Thousands of hours of computing time are 
spent identifying and removing RFI from the datasets, before one can search for SETI 
signals. Since its inception, the SETI@home project has detected over 4.2 billion potential 
signals. While essentially al l these potential sig nals are RFI, it is possible that therein lies a 
true extraterrestrial signal ([Korpela et al.ll2010l ). VLBI offers several ways of discriminating 
RFI signals from SETI signals not applicable to single dish or even short baseline radio 
interferometers. 

First, VLBI uses N telescopes separated by hundreds to thousands of kilometres giving 
N(N — l)/2 baselines, once the data from independent pairs of telescopes are correlated 
together. This aspect of VLBI is important, as RFI which is not present at both telescopes 
involved in a baseline does not correlate. An astronomical or SETI compact radio source 
within the interferometer field of view will be detected on all baselines. This is also true for 
RFI if it is detectable by both telescopes in a baseline. However, even correlated RFI can 
be easily identified with VLBI, as it produces distinct phase variations in both frequency 
and time in the correlated output as discussed below. 



3.3. VLBI Phase Variations 

The correlated signals from an interferometer are complex quantities (known as 
visibilities), with amplitude and phase as the fundamental observables. The amplitudes 
are proportional to the correlated power of the source, while the phase, <fi represents the 
residual time delay that results from the wavefront of the radio emission arriving at one 
antenna before the other. At the centre of the field of view (also known as the phase 
centre), the delay (and hence the phase) is set to zero, via application of a a priori time 
delay, r, at the correlator. This time delay accounts for many different effects, but primarily 
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simple geometry. The phase, of a radio point source offset from the phase centre can then 
be related to i ts pos ition relative from the phase centre, and is given by equation 1 from 



Morgan et al 



(2011). 



2nv 



(lu + mv) 



(1) 



where v is the observing frequency, I and m are the coordinates of the source in the 
plane of the sky, relative to the phase centre, and u and v are coordinates which describe 
the instantaneous geome try of the projected ba seline , with respec t to th e phase centre. 



The reader is referred to 



Thompson et al 



Tavlor et all (119991 ) for details. For a 



(120011 ) and 

finite observational bandwidth, and for the duration of the observation, the radio source in 
equation [TJ will display a phase change with frequency and time given by, 



d(p 2tt 

— = — (tit + mv), 
au c 



(2) 



dd> 2ixv ,,du dv. 

-r = (l-r + 

dt c dt dt 



(3) 



By measuring the phase change with time and frequency for a given source, its location 
with respect to the phase centre can be determined using the above equations. Additionally, 
if the source's position is known, maximum values for both d<fi/dv and dcj)/dt can be 
determined, using maximum values of u or v and du/dt or dv/dt. 



3.3.1. Determining the Phase Variations with Fourier Transforms 



An estimate of d<fi/dt or d<fi/dv for a radio source can be determined by the discrete 
Fourier transform (DFT) of the complex visibilities. For narrow-band (persistent or 
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intermittent in time) and broad-band signals, DFTs can be applied to the complex 
visibilities as a function of time to determine d(j)/dt and frequency to determine d<fi/di/, 
respectively. 

This method is illustrated in Fig. [TJ The upper panel shows the phase of a narrow-band 
test signal as a function of time, with a sample time of 5 minutes per point over a 40 minute 
observation. The lower panel displays the corresponding discrete power spectrum obtained 
by applying DFTs to the complex visibilities. 

A signal with little or no phase change in time or frequency (implying a source close 
to the phase centre) will appear in bin zero. However, a signal with high phase change 
with time (as the test signal in Fig. [TJ) is indicative of a source far from the phase centre, 
probably RFI, and will appear in a different bin. The height of the DFT amplitude of bin 
zero compared with the others gives an indication of the signal to noise of any signal which 
comes from the phase centre. 

There are, however, ambiguities in signals with d<fi/dt greater than the Nyquist 
frequency. The component of the test signal located in the —3 bin can be generalised as 
dcj)/dt = —135 ± 180iV o per data point, where iV is a a wide range of positive and negative 
integers. The d(f>/dt of the Nyquist bin, due to ambiguities is ±180° per data pointy. 



2 The reader is refereed to iBracewelll ( 120001 ) for further details on discrete Fourier trans- 
forms 



- 10 - 



4. Searching for Signatures of Extra- Terrestrial Intelligence from Gliese 581 



4.1. Observations 



G1581 was observed on 2007 June 19, for 8 hours in the 20 cm wavelength band 
by three stations of the Long Baseline Array (LBA): Mopra (Mp), Parkes (Pa) and the 
Australia Telescope Compact Array (At) in phased array mode. The observation was 
conducted in three frequency pairs of 64 MHz bandwidth centred on: 1262,1312; 1362,1412; 
1462,1512 MHz, in dual-polarisation mode. The uv coverage is shown in Figure [2j This 
configuration gives a maximum resolution of 127 mas. The observations of G1581 were 
phase referenced to the ICRF2 source 1510-089 (J1512-0905, RA = 15 fe 12 m 50.532925 s ; De c 



-09°05'59. 82961", J2000), a calibrator located 2.14 degrees from G1581 flFev et al. 



200J). 



1510-089, shows an unresolved core with extended structure to the North 



at 2.3 GHz (IFev fc Chariot 



19971), and 327 MHz (IRampadarath et al.ll2009l ) on the longest 



last of the core 



baselines of the Very Long Baseline Array (VLB A). The longest baseline of our observation 
(At-Pa) is a factor of 9 less than the longest VLBA baseline, hence, 1510-089 is expected to 
be unresolved in our observations. 

A phase-referencing cycle of 5 minutes (2 minutes on calibrator, 3 minutes on G1581) 
was employed. A single observation block consisted of 8x5 minute phase-referencing cycles. 
The entire observation comprised 12x40 minute observational blocks switching between the 
three frequency pairs (table 1). 



20071 . 



2011) 



The data were correlated using the DiFX software correlator (iDeller et al 
with a frequency resolution of 1.95 kHz per spectral channel, and converted to the standard 
FITS-IDI format. The lowest and highest 5 MHz (2500 channels) were discarded from each 
band due to band edge drop-off in sensitivity. The data files were then sub-divided into five 
datasets of 4000 channels each and one dataset of 1268 channels at each centre frequency. 
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4.2. G1581 Coordinates 

Due to its proximity to Earth, the exact coordinates of G1581 at the time of observation 
must be calculated taking into account parallax and proper motion. An accurate position 
was computed via a python script (James Miller- Jones, ICRAR, private communication), 
using a reference position of G1581 (RA and Dec) and the barycentric coordinates of the 
Earth (in AU) at the time of the observation (obtained from NASA's Jet Propulsion 
Laboratory HORIZONS on-line solar system data and ephemeris computation servicqj). 



The reference position for G1581 was taken from the Hipparcos catalog ( IPerryman et al. 



19971 ). G1581 was observed on 1991 March 1, by the Hipparcos satellite, and a position 
of RA(J2000)=15 ?i 19 m 26.825 s ; DEC(J2000)=-07°43'20.210", was recorded with positional 
errors: e#A=l-99 mas and e£>,Ec=l-34 mas. Additionally, a parallax of 159.52 ± 2.27 mas 
was measured by Hipparcos. A proper motion of RA = —1227.67 ± 3.54 mas/yr; DEC 
= —97.78 ± 2.43 mas/yr, was also taken from the Hipparcos catalog. The position 
calculated for G1581 at the VLBI observation date, is RA(J2000)=15 h 19 m 26.206 s ; 
DEC(J2000)=— 07°43'17.555". This calculation has a positional uncertainty of 53 mas, 
largely due to the cumulative effect of the uncertainty in the proper motion measured by 
the Hipparcos 1991 observation. 



4.3. Data Reduction and Calibration 



The data reduction and calibration was performed using the National Radio 
Astronomical Observatory (NRAO) software package AIP9tj. Prior to calibration and data 



3 http: / / ssd.jpl.nasa.gov/horizons.cgi 

4 The Astronomical Image Processing System (AIPS) was developed and is maintained by 
the National Radio Astronomy Observatory, which is operated by Associated Universities, 
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analysis, the phase centre of the calibrator was shifted, via the AIPS task UVFIX, to match 
the calculated positional difference. Flagging and phase-only calibration (fringe-fitting 
and self-calibration) were carried out on the phase centre shifted calibrator. The phase 
solutions were transferred and applied to the G1581 datasets. After calibration, further 



analysis was carried out using the ParselTongue scripting language (jKettenis et al. 



2006|) 



and the numerical python package, numpyo. The flagging, phase calibrations and transfer 
of solutions were implemented separately for the frequency pairs listed in Table 1. 



4.4. Expected Upper limits to dcfi/dv and d<p/dt for the G1581 system 

The uncertainty in the position of G1581 (53 mas) and the semi-major axes of the 
planets (G1581d = 40 mas and G1581g = 23 mas) produces a quadratic error of 66 mas and 
57 mas on the positions of G1581d and G1581g, respectively. The positional uncertainty 
places upper limits on the dip/dv and d<p/dt expected for a SETI signal from G1581. From 
equations [2] & [3] the upper limits are: (d(j)/ 'du) max = (1.2 x 10~ 10 )° per spectral channel 
and (d(p/dt) max = (2 x 10~ 4 )°per second (or (6 x 10~ 2 )°per data point, after averaging 
each scan). 

Inc., under cooperative agreement with the National Science Foundation. 
5 http:/ /numpy. scipy.org/ 
6 1 spectral channel = 1.95 kHz 
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4.5. Analysis 

4-5.1. Search for Candidate SETI Signals 

After correlation and calibration our dataset consisted of 8 hours of data for the three 
baselines, divided into 12x40 minute observational blocks. Each 40 minute observation 
block was divided into three minute scans of G1581, consisting of two second time- averaged 
complex visibilities and 32768 spectral channels. Before any further analysis was performed, 
each scan was vector-averaged along the time axis (for three minutes) leaving the 
observational blocks as arrays of data, with 32768 channels and 8 time integrations. 

Since the star G1581 is not a radio source (see Sect. [2]), the dataset was not expected 
to contain any continuum emission. The dataset therefore consisted almost entirely of 
Gaussian noise at the thermal noise level with a few narrow spectral lines, possibly transient 
RFI. Any SETI signals present are expected to be weak compared to RFI. For data 
containing only weak signals (signal-to-nose ratio, SNR <C 1), the mean of the amplitudes 



of th e complex visibilities can be used as a proxy for the thermal noise ( iThompson et al. 



2001 



, Sect. 9.3). However, the mean is biased by RFI. To avoid this bias, the median of 
the amplitudes of the complex visibilities was used in place of the mean as a proxy for the 
thermal noise. 



Broad-band signals 

Each vector-averaged scan was searched along the channel axis for visibilities with 
sufficient amplitudes to have a signal-to-noise ratio (SNR) > 5. Detected signals were then 
checked to see if they occur on all baselines. Signals found on all baselines were plotted 
as a function of channel and searched for instances of five or more adjacent channels with 
amplitude and phases displaying characteristics of broad-band signals. There were 22 
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groups of signals satisfying these requirements, selected as broad-band SETI candidates. 



Narrow-band signals 

The observational blocks, consisting of 32768 channels and 8 time integrations, were 
vector-averaged along the time axis, leaving 1-dimensional datasets of 32768 channels. The 
data were then searched for visibilities with sufficient amplitudes to have an SNR > 5. 
Detected visibilities were then checked to see if they occur on all baselines. There were 200 
signals satisfying these requirements, selected as narrow-band SETI candidates. 



Frequency Distribution of the Candidate Signals 



Figure |3] shows the frequency distribution of the SETI candidate signals (1525- 
1534 MHz), with the 200 narrow-band (grey) and the 22 broad-band (black) candidate 
signals. Australian Space to Earth geostationary satellites, Optus' Mobilesat and the 



INMARSAT 



3 based systems, marketed i n Austra 



frequencies in the range 1525-1559 MHz (jSymons 



ia by Telstra are known to operate at 



19981 ). The close match between the 



transmission frequencies of the satellites and our candidate signals suggests that many or 
all of our detected signals may be RFI coming from these satellites. 



7 The International Mobile Satellite Organisation (INMARSAT) is an international treaty 
organisation with sixty seven member countries including Australia. INMARSAT operates 
a GSO satellite system for global, mobile communications on land, sea and in the air 



- 15 - 



4-5.2. Phase Variation of broad-band Candidate SETI Signals 

Recall from Sect. I3.3[ dtfi/dv is dependant on the location of the signal with respect 
to the phase centre. To obtain d^jdv for the 22 broad-band SETI candidates, DFTs were 
applied to the complex visibilities as a function of frequency, as described in Sect. 13.3.11 
The maximum expected d<f)/dv for a radio signal from the planets is (1.2 x 10~ 10 )°per 
spectral channel (Sect. 14. 4p . This phase variation is effectively zero at the level of frequency 
sampling in the observations. 

Figure H] shows the d^jdv for the 22 broad-band SETI candidates. The horizontal 
dashed lines indicate the d(f)/dv range expected for a radio signal from the G1581 system. 
Any signal with a d<p/dv range that crosses the horizontal dashed lines on all three 
baselines are considered as possibly originating from the G1581 system, and merits further 
investigation. None of the 22 broad-band signals has satisfied this criteria, and are thus 
rejected as originating from the G1581 planets. 

A large percentage of the signals are clustered around dtfi/dv = 45° per spectral 
channel on the shortest baseline (At-Mp) in Fig. HI Using equation |2l and assuming a 
North-South baseline of projected length, v = 0.5 MA (see Fig. [2]), this phase variation 
implies a minimum distance of ~ 10° from the phase centre (c.f. declination of G1581 ~ -7° , 
geostationary satellite declination ~ +5°). This effect is less clear on the longer baselines 
since the rapidly rotating phase of the complex visibilities along both the frequency and 
time axes leads to averaging losses (smearing). 

4-5.3. Phase Variation of Narrow-Band Candidate SETI Signals 

The discrete power spectra of the 200 narrow-band SETI candidate signals were 
obtained by application of DFTs to the complex visibilities as a function of time. As 
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discussed in Sect. I4.4[ the maximum d(f)/dt expected from a radio signal originating from 
either of the planets is 2 x 10~ 4 degrees per second. Such a slow varying signal is expected 
to behave as a signal from the phase centre, where most of its power appears in bin zero of 
the DFT power spectrum on the three baselines. 

To determine whether any of the signals are from the G1581 system, classical detection 
theory was employed to test the hypothesis that the signal is located at the phase centre. 
Hypothesis testing compares the likelihood that any given signal is consistent with a source 
located at the phase centre, to the likelihood that it is consistent with a source located away 
from the phase centre, or from RFI. This i s the class ical Neyman-Pearson (NP) approach 



to hypothesis testing, or signal detection (IKay 



1998 



chap. 3). The classical detection 



method has the advantages of providing robust estimation of errors and an efficient method 
of discriminating false positives from true signals. 



Hypothesis Testing and Monte-Carlo Simulations 

For Gaussian-distributed noise in the complex visibilities, the amplitudes of the 
DFT bins are Rician-distributed. The two hypotheses are as follows: (1) The alternative 
hypothesis, r K\\ signal-present - a source located at the phase centre has most of the 
power concentrated in bin zero, yielding a Rice distribution of amplitudes with parameters 
(amplitude and variance) determined by the data SNR and noise level; (2) The null 
hypothesis, Hq] signal-absent - a source located away from the phase centre, yielding power 
in a bin other than bin zero, with the same Rician parameters. The ratio of the likelihoods 
for bin zero to the next highest bin is the test statistic. 

The test statistic is compared to a threshold to determine our level of confidence that 
the signal is real. To determine the threshold, we simulate 10 6 signal-present and 10 6 
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signal-absent sources, with the non-phase centre signals distributed randomly in the field, 
and for SNRs that reflect that found in our dataset. In addition, given that we expect the 
primary signals in the data to be RFI, rather than true signals, we also simulate zero-mean 
Gaussian noise with random RFI spikes located at a single time point in the visibilities. 
Figure |5] plots the log-scale distribution of the test statistics for both the signal-present 
(upper plot) and the signal-absent case (lower plot). T he thresho lds were chosen so as to 
minimise the chance of false-positives, i.e Type I error (jKaylll998l . chap. 3). From Figure |5] 
the thresholds ranges from the maximum log of the test statistic for the signal-absent case 
(= 2), and the minimum log of the test statistic for the signal-present case (= 35). 

The method was then applied to the 200 narrow-band signals, where the overall test 
statistic, T(x) was taken as the product of the test statistics of the individual baselines. 
Ten narrow-band signals were found to have ln[T(x)] > 2, of which only two are beyond 
ln[T(s)] = 35 (Table 2). The data are expected to be consistent with no signal present 
at the phase centre (i.e. consistent with the null hypothesis). The significance level of 
rejecting the null hypothesis (the p-value) for the 200 signals were found by integrating the 
probability distribution of the signal-absent test statistic, obtained from the simulations. 
A small p-value (p < 0.01) corresponds to a signal that may not be consistent with the 
null hypothesis, and warrants further investigation. Figure E] shows the p- values of the 200 
signals. Most of the 200 signals were found to have a p-value greater than 0.01, with ten 
signals having a p-value of < 0.01. 

Eight signals lie between both distributions (2 < ln[T(x)] < 35). These signals are 
strictly inconsistent with both hypotheses, and reflect the limitations of our simple analysis; 
namely the simulation of data with either pure signals at the phase centre, or a single RFI 
peak in the visibility, and the need for the overall test statistic to have a value in excess of 
the threshold, rather than those for the individual baselines. 
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Table 2 lists ln[T(a;)] along with the test statistic for the individual baselines for all 10 
events with ln[T(x)] > 2. For all but one of the 10 events, ln[T(x)] < 2 if the shortest 
baseline is excluded. For all 10 events the test statistic is far greater for the shorter baseline 
than for the other baselines. This is not consistent with a source at the phase centre. It 
is, however, consistent with RFI distant from the phase centre. We can therefore easily 
exclude these signals as being from the GL581 system. 



5. Discussion 

We present the first VLBI SETI experiment. This pilot study demonstrates that 
the fundamental properties of VLBI makes it an ideal technique for targeted SETI. With 
multiple baselines, milliarcsecond resolution, and millijansky sensitivities, RFI can be easily 
identified and rejected using automated techniques, while being sensitive to very weak 
signals. 

Our search method targeted signals that are five times the baseline sensitivity, which 
were then cross-referenced over three baselines, to identify corresponding signals. Hence, 
the amplitude sensitivity of our array was limited to our least sensitive baseline over an 8x3 
minute observing interval, for narrow-band signals (At-Mp; la = 0.31 mJy) and a single 3 
minute scan for broad-band signals (a factor of \/8 less sensitive jf. 

Our experiment finds no evidence for radio signals > 1.55 mJy (= 5cr) within the 
frequency range of 1230 - 1544 MHz from the region of G1581. We can therefore place an 
upper limit of 7 MWHz -1 on the power output of any isotropic emitter located in this 
planetary system, within this frequency range. On the other hand the power of the Arecibo 



8 The sensitivity limits were derived using the System Equivalent Flux Densities (SEFD) 
for the individual telescopes at 1400 MHz; At=68 Jy, Pa=40 Jy, Mp=340 Jy. 
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planetary radar transmitter is ~ 1 MW (IBlack 



2002 



Tarter 



20041 ). Assuming an aperture 



efficiency of 50%, at A = 2 1 cm the directive power output is ~6 TW. Over a bandwidth 



of 5.4 Hz (Black et al. 



20011 ). this emission would have the flux density of a ~200 Jy radio 



source (or 650 mJy if integrated over our channel width of 1.93 kHz) at a distance of 20 ly. 



Using current VLBI instruments for SETI 

Since this was a proof-of-concept observation, two differences between this experiment 
and a typical VLBI observation should be noted. Firstly the baseline lengths in this 
experiment were short by VLBI standards. Baselines one or even two orders of magnitude 
longer would significantly reduce the amount of correlated RFI. Secondly, most VLBI 
observations are carried out using more than three telescopes. Increasing the number 
of baselines to four or more constrains both the phase and amplitude closure relational, 
drastically increasing the quality of calibration. This would also provide the opportunity 
for image plane searching. However, it would likely be computationally expensive as a 
technique for searching over very wide frequency ranges. Relaxing the single baseline five 
times median limit could enhance the sensitivity beyond that of the weakest baseline. 

We should also note that the width of the DFT bin is not a fundamental limit on 
the determination of the phase variation. A least-squares fit would allow the delay and 
delay rate to be determined more accurately. Essentially, the error on these parameters 



is limited by the signal t o noise. Certainly 



sub-milliarcsecond level ( iDeller et al. 



;he LBA has been used for astrometry at the 



2009b 



a). Observations of spacecraft have proven 



9 The sum of the phases on any triangle of three baselines, and the combined amplitude 
ratio of four baselines c ontains information onl y on the true visibility of the source itself, all 
other errors cancel out ( [Thompson et al.ll200ll . chap. 9). 
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VLBI astrometry techniques for very weak narrow-band sources (lAvruch et al.ll2006l ) and 
suggest that VLBI might even allow the detection of the orbital motion of SETI signals. 

While the techniques discussed in this paper are ideal for targeted SETI, we also note 
the possibilit y of carrying out SETI commensally with normal VLBI operations. SETI 



target stars (ITurnbull fc Tarter 



20031 ) would be expected to be present in the field of view 



of the antennas for V LBI observations reasonab ly regularly. Using the techniques discussed 



in 



Deller et al 



(120111 ) and 



Morgan et al 



(1201 if ) it would be possible to produce visibility 
datasets for any target star in the primary beam. Additionally, direct analysis of the voltages 
(the baseband data) although computationally expensive, could also be implemented at the 
correlator in an efficient manner to search for signals from extraterrestrial intelligence. A 
possible model for this is the VFASTR transient search cur rently running at the VLBA 



operations centre (jWayth et al 



2011 



Thompson et al 



2011) 



SETI and the Square Kilometre Array (SKA) 



A sensitive and advanced radio telescope, the Square Kilometre Array (SKA) is 
currently being planned for the coming decade^. 



At 1400 MHz , the long baselines of t he SKA will provide an astrometric precision of 



3 yU-arcseconds (IFomalont fc R eid 



sub-/xJy sensitivities ( ICarilli &: Rawlings 



200 



. In addition, the SKA is expected to achieve 



20041 ) . At ~20 ly the SKA would probe an isotropic 



luminosit y of a few kW Hz , far below the power output of radars and telecommunication 



satellites rtTarter 



20041 ) . The fields of view will be ~1 sq uare degree at the maxim um 



baselines and ~30 square degrees for the shorter baselines (ICarilli &: Rawlings 



20041 ). These 



capabilities will make the SKA the fastest, most sensitive and widest field of view survey 



10 



http : //www . skatelescope . org/about/pro j ect/ 



instrument in astronomy. However, for the SKA to be used for SETI, suita ble targets are 
required. Such targets ar e currently being provided by the Kepler Mission (IBorucki et al. 
201ll ; Batalha et al-lhoij ). 



6. Conclusion 

The discoveries of potential habitable planets are no doubt fuelling a renewed 
excitement in SETI. With current VLBI systems we are able to search these candidates for 
signals from extra-terrestrial civilizations over a wide range of frequencies, and with high 
sensitivity. 

In addition, VLBI is undergoing a transition from traditional recording systems towards 
real-time correlation via network connections (so called e-VLBI). This has led to VLBI 
being used for rapid follow-up of transient sources. These developments, the wide range of 
frequencies covered by VLBI instruments, and the huge added value that VLBI techniques 
bring to a targeted SETI search make VLBI the natural choice for follow-up observations 
to confirm a potential SETI signal. However, more accurate positions are needed for target 
stars for future experiments with higher angular resolution. Such p ositions will be provid ed 



in th e near future by the high-quality optical astrometry of Gaia ( IBastian Schilbach 



1996|). 
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the Australia Telescope which is funded by the Commonwealth of Australia for operation 
as a National Facility managed by CSIRO. SJT is a Western Australian Premier's Research 
Fellow, funded by the state government of Western Australia. The Centre for All-sky 
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UT 


T 1 T "\ • /TV ITT \ 

rrequency Pairs (MHz) 


07:40:00-08:30:00 


1262,1312 


08:30:00-09:10:00 


1362,1412 


09:10:00-09:50:00 


1462,1512 


09:50:00-10:30:00 


1262,1312 


10:30:00-11:10:00 


1362,1412 


11:10:00-11:50:00 


1462,1512 


11:50:00-12:30:00 


1262,1312 


Iz.oU.UU-lo. 1U.UU 


1 QfiO 1/110 


13:10:00-13:50:00 


1462,1512 


13:50:00-14:30:00 


1262,1312 


14:30:00-15:10:00 


1330,1380 


15:10:00-15:30:00 


1462,1512 



Table 1: The observation was sub-divided into 12 separate sessions, that alternated between 
3 frequency pairs. This table lists the observation times and the corresponding frequency 
pairs. 
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Frequency 

1 J 


UT 


L \ / J 


Baseline Specific log of the test statistic 


(MHz) 






At-Mp 


At-Pa 


Pa-Mp 


1531.69080115 


13:10:00-13:50:00 


2.39 


1.92 


-0.30 


0.77 


1533.6297203 


09:10:00-09:50:00 


3.11 


6.12 


-1.36 


-1.65 


1533.80206867 


13:10:00-13:50:00 


3.50 


2.19 


0.66 


0.66 


1527.85017443 


09:10:00-09:50:00 


4.21 


4.89 


1.40 


-2.07 


1527.50351919 


11:10:00-11:50:00 


6.51 


16.91 


-4.95 


-5.45 


1527.56227431 


11:10:00-11:50:00 


7.25 


11.04 


-0.27 


-3.52 


1525.18856723 


09:10:00-09:50:00 


14.84 


30.53 


-15.32 


-0.37 


1527.50156068 


09:10:00-09:50:00 


26.19 


17.86 


2.87 


5.47 


1532.03158088 


09:10:00-09:50:00 


86.19 


88.270 


0.03 


-2.10 


1527.56423282 


09:10:00-09:50:00 


242.77 


1782.83 


-1305.76 


-234.30 



Table 2: The ten signals with ln[T(x)] > 2. Column 1: The centre frequency; Column 2: 
The observed (UT) time period; Column 3: The sum of the log of the test statistics of the 
individual baselines, which is taken as the final test statistic; Columns The log of the 
test statistics of the individual baselines. 
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DFT bins 



Fig. 1. — The phase variation with time (upper panel), and the corresponding discrete 
Fourier transform (DFT) for a narrow-band test signal. The DFT bins provide an estimate 
of the d(f)/dt of the signal. The Nyquist bin at bin ±4, represents the maximum possible 
d(j)/dt, 180° per data point. Each data point has a sampling time of 5 minutes. The test 
signal is found to have a d<f>/dt of —135° per point, equivalent to 3 turns over the time range. 
However, due to ambiguities d<fi/dt of this signal can be generalised to — 135±180iV o per data 
point, where N is a wide range of positive and negative integers. The same principle applies 
in determining d(f)/du for broad-band signals, where the DFT is applied to the complex 
visibilities as a function of frequency. 
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Fig. 2. — uv coverage of the LBA observation of the G1581 system, in units of mega wave- 
lengths. The x-axis plots the u coordinates, while the v coordinates is on the y-axis. 
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Frequency/MHz 

Fig. 3. — Distribution of the SETI candidate signals, 200 narrow-band emission (grey) and 22 
broad-band events (black). The SETI candidate signals fall within the operating frequency 
range (1525-1559 MHz) of the geostationary Australian Space-Earth transmission satellites, 
Inmarsat-Optus. 
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Fig. 4. — The dip/dv range and bin number for the 22 SETI broad-band candidate signals, 
for the three baselines is plotted on the y-axis {upper panel: At-Mp; centre panel: At-Pa; 
lower panel: Mp-Pa). The error bars give the range of possible d(p/di> (or bin widths) for 
each signal. The x-axis gives the event number (1-22). The dashed line marks the maximum 
dtp/dv expected for a radio signal from G1581, effectively zero at this resolution. 
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Fig. 5. — Log-scale distribution of test statistics, ln[T(x)] obtained from Monte-Carlo simu- 
lations for both I. Signal-absent (upper plot), and 2. Signal-present (lower plot). Note the 
very different ranges plotted on the x-axis. 
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Fig. 6. — Distribution of p-values for the 200 narrow-band signals. The p-values represent 
the error in rejecting the null hypothesis, , "H (signal-absent). The line plots the p = 0.01 
significance value. Signals at or below this line can be considered inconsistent with the signal- 
absent hypothesis, and require further analysis. The signals at p=l all have test statistics 
lower than the minimum test statistic of the simulated signal-absent distribution (Fig. [5]), 
and are assigned a p- value of 1. 



